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Abstract 

The polarization of Polyakov type loops is responsible for the dif- 
ference between quenched and unquenched finite size effects on the 
QCD mass spectrum. With a numerical simulation, using different 
sea quarks boundary conditions, we show that we can align the spa- 
tial Polyakov loops in a predefined direction. Starting from these 
results, we propose a procedure to partially remove the Polyakov type 
contributions in the meson propagators. 



1 



The finite extent of the lattice is an important source of systematic errors 
in lattice QCD calculations. Theoretical and numerical analysis |2[ faced 
recently the problem of finite size effects in full lattice QCD. The conclusions 
of these analysis are that the behaviour of the hadronic masses as a function 
of the lattice size L follows a power law uil = + cL~ u , that only asymp- 
totically goes to the exponential decay predicted by the effect of virtual pions 
emitted from a point-like hadron || . 

The magnitude of finite-size effects is much smaller for quenched QCD 
than for full QCD, especially below La ~ lfm: v =1-2 for the quenched case 
and v =2-3 for full QCD. 

The reason of the difference can be understood using, for example, the 
hopping parameter expansion on the valence quark mass that leads to the 
following relation for the meson propagator ||: 

E ^\W{C)) + £ k l v ^a val (P(C)} (1) 
c c 

where the sums extend over all possible closed paths (C) of length 1(C), the 
W(C) are the Wilson loops completely contained into the lattice while P(C) 
are valence quark loops wrapping around the lattice in the spatial directions 
(Polyakov-type) and (•) denotes field averages; the a va i represents the spatial 
boundary conditions on the valence quarks: a va \ = +1 for the periodic and 
&vai = — 1 for the antiperiodic cases. 

The averaged Polyakov loop (P) is different from zero in full QCD, while 
it is zero in the confined phase of quenched QCD. This means that the second 
term in eq.(l), due to the loops which go around the lattice, gives different 
finite size effects between the quenched and the full QCD value of mesonic 
correlations. 

To obtain finite size effects comparable with those of the quenched case 
we have to remove the Polyakov loop contributions. 

In this paper we want to demonstrate that we can align the Polyakov 
type loops in a predefined direction using sea quarks boundary conditions. 
Starting from these results we propose a procedure to remove the second 
term of eq.(l). 

On a finite lattice with periodic boundary conditions on the gauge fields 
there is a symmetry in the pure gauge action consisting in multiplying all 
links in the /i direction at a constant x^-plane by the complex number z^, 
that belongs to the centre of the gauge group (Z 3 )j4|: 
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z k = e 12 ^ 3 v k = (0,1,2). (2) 

Under it the Polyakov loops in the \i direction are not invariant, but they 
transform as 

P - z k P (3) 

In full QCD we consider both the gauge and the fermionic action. In the 
fermionic action 

Swuson = -k^ii'i^i 1 -li^U^x) ifj(x + ii) 

X,fl 

+1>(x)(l + ln)U}(x -n)ij}{x-n)) ^ 

x 

the kinetic part is not invariant under Z%, since on the boundary it is not 
possible to find a gauge transformation of the fermionic fields that cancels 
out the centre transformation on the gauge fields. Thus the symmetry, that 
in the quenched confined case guarantees that (P) = 0, is explicitly broken 
by the kinetic part of the fermionic action. Because the non-invariant term is 
proportional to k the amplitude of this violation is more important for light 
sea quarks. 

It is possible to summarize what happens on the lattice boundary. Making 
a double expansion in the full QCD, firstly with a strong coupling expansion 
on (3 and secondly with a hopping parameter expansion on the sea quark 
masses, we obtain the 3-d Potts Model with magnetic field. The introduction 
of the fermionic action in QCD is equivalent to turn on a magnetic field 
h which breaks the Z 3 symmetry. To study the situation on the lattice 
boundary we introduce a simple model of a single spin II that can take the 
three possible values: 

n = 1, III = e l2 " /3 , n 2 = e- l2n/3 (5) 
with a Hamiltonian 

H = hU + /i t n t (6) 

which for h ^ is not Z 3 invariant and it is composed by terms like those 
which break the Z% simmetry in the Potts Model with magnetic field. 
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The II argument is related to the phase of Polyakov loop and the value 
of h with the sea fermionic boundary conditions. 

We summarize the interesting h choices in fig. la-le. We see that with 
h = +\h\ (periodic boundary conditions on sea quarks) there are two pre- 
ferred states, while with h = —\h\ (antiperiodic boundary conditions on sea 
quarks) only the II = 1 state is selected. This circumstance suggests a way 
to select the other two Z 3 states, as we can see from fig. Id and le. 

Thus we expect that with periodic boundary conditions on sea quarks 
both the e l2?r//3 and the e~ l2n ^ 3 phases of the Polyakov loops are preferred 
by the system. This means that for each gauge configuration there will be 
domains in which they will point towards either one of the two directions 
e i27r/3 anc j e -«2?r/3_ suc j 1 structures exist with the antiperiodic boundary 
conditions and moreover Polyakov loops are likely to point towards 1 in the 
Z 3 space. 

Moreover we can align the Polyakov loop in the e l2n ^ 3 (or e _l27r//3 ) in the 
Z3 space if we choose — e~ t2n ^ 3 (or — e t27r / 3 ) boundary conditions on the sea 
quarks. 

To check the foregoing suggestion we performed on APE100 a full QCD 
simulation I| with 2 flavors Wilson fermions at j3 — 5.3 on a 8 3 x 32 lat- 
tice with k sea = 0.1670. We performed two different runs, one with peri- 
odic boundary conditions on the sea quarks and the other with antiperiodic 
boundary conditions. We collect in both cases 600 thermalization trajectories 
plus other 1200. On the latter we perform, every 5 trajectories, a measure- 
ment of the spatial Polyakov loops. We use the smearing procedure for the 
measurement of Polyakov loops || for 10 values of smearing. 

The situation for the phases of the spatial Polyakov loops are reported 
in figures 2 and 3. With antiperiodic boundary condition, fig.2, we obtain 
that the phase is close to zero. Otherwise with periodic boundary conditions, 
fig. 3, we obtain that the phases are spread in regions near e l27T ^ 3 and e~ l2n ^ 3 . 

We also impose the boundary condition — e~ l2n ^ 3 and — e t2n ^ 3 on the sea 
quarks in a quick simulation on a 4 3 x 6 lattice with (3 = 3.0 and k sea = 0.1670 
and we verified that Polyakov loops point towards the e* 2?r//3 and e~ l27T ^ 3 
correspondingly, see fig. 4. 

From eq.(l) (see ref. f2j), we do not expect large differences between the 
meson masses calculated either with periodic boundary conditions on sea 
and valence quarks or antiperiodic boundary conditions on sea and valence 
quarks, because the sign (not the amplitude) of the second term is the same 
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in both cases. But due to our results on the polarization of the Polyakov 
loops we conclude that the best boundary conditions are the antiperiodic, 
because their Polyakov loops have lower dispersion than with the periodic 
ones. 

Moreover the previous analysis gives us a procedure to partially eliminate 
the Polyakov loops that contribute to eq.(l). In fact we can choose the an- 
tiperiodic boundary conditions on the sea quarks. In this way all the gauge 
configurations have Polyakov loops that point towards the 1 Z(3) state. Then 
we invert three times the fermionic valence operator: with antiperiodic, with 
_ e i27r/3 anc j _ e -«27r/3^ rj^g avera g e G f the meson propagators obtained 

with the three inversions has a factor a va i = in the second term of eq.(l). 
We note that the required CPU time for valence fermionic inversion is negli- 
gible compared to the CPU time needed to obtain the full QCD configuration 
and that in the full QCD simulations the CPU time scales as Volume 5 ^ 4 . 

The numerical simulations of this work have been done with configura- 
tions obtained using 2 months of CPU time of a 128 nodes APE100 machine. 
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FIGURE CAPTIONS 



Figure 1. The energy levels of the model of eq.(6) for different values of h. 
The notation for the states of II is that of eq.(5). 

Figure 2. Histogram of the phase of the x, y and z components of the 
Polyakov loop for antiperiodic boundary conditions on the sea quarks. Data 
are from trajectory 440 to trajectory 1800. The value of smearing is 10. The 
lattice is 8 3 x 32. 

Figure 3. The same of Figure 2. for periodic boundary conditions on the 
sea quarks. 

Figure 4. Behavior of the phase of the average of the three spatial Polyakov 
loops as a function of the HCMA trajectories. Data are from trajectory 205 
to trajectory 600. The value of smearing is 0. The lattice is 4 3 x 6. In 
fig. 4a we report the case of periodic boundary conditions on the sea quarks; 
in fig. 4b we report the case of antiperiodic boundary conditions on the sea 
quarks; fig. 4c is obtained from the boundary condition — e~ l27r / 3 and fig.4d 
from the boundary condition — e t2n ^ 3 . 
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